Constraint on the growth factor of the cosmic structure from the damping of the 

baryon acoustic oscillation signature 
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We determine a constraint on the growth factor by measuring the damping of the baryon acoustic 
oscillations in the matter power spectrum using the Sloan Digital Sky Survey luminous red galaxy 
sample. The damping of the BAO is detected at the one sigma level. We obtain o&Di{z = 0.3) = 
0.42^Q 2g at the la statistical level, where erg is the root mean square overdensity in a sphere of 
radius 8/i _1 Mpc and D\{z) is the growth factor at redshift z. The above result assumes that other 
parameters are fixed and the cosmology is taken to be a spatially flat cold dark matter universe 
with the cosmological constant. 

PACS numbers: 98.80.-k,95.35.+d,95.36.+x 



I. INTRODUCTION 

The baryon acoustic oscillations (BAO) are the sound 
oscillations of the primeval baryon-photon fluid prior to 
the recombination epoch. The BAO signature imprinted 
in the matter power spectrum is very useful for the study 
of the dark energy, hypothetically introduced to explain 
the accelerated expansion of the universe [J, 0| , because 
the characteristic scale of the BAO plays a role of a stan- 
dard ruler in the universe [1, 0, @ . The BAO signature 
in the galaxy clustering is clearly detected [H, 0], and 
the constraints on the equation of state parameter of the 
dark energy are demonstrated [j| H, EH EH • Future large 
surveys for the precise measurement of the BAO are in 
progress or planned, providing us an important tool to 
explore the origin of the accelerated expansion of the uni- 
verse. 

Besides the measurement of the BAO, those future sur- 
veys provide us other important information. For exam- 
ple, the redshift-space distortions will be measured pre- 
cisely at the same time. The redshift-space distortions 
reflect the velocity of galaxies in the direction of the line 
of sight. Especially, the linear redshift-space distortion 
due to the Kaiser effect comes from the linear velocity 
field of matter perturbations, whose measurement will 
give us a chance to test t he g ravity theory on the cosmo- 
logical scales [H, EH EH Il5f This is important because 
several modified gravity models are proposed to explain 
the accelerated expansion of the universe as an alterna- 
tive to the dark energy. Thus the measurement of the 
redshift-space distortion, which is useful to constrain the 
growth rate and the growth factor of the cosmic struc- 
ture, is also important. 

In Ref . E3] , it is pointed out that the precise mea- 
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surement of the BAO may also be useful to constrain the 
growth factor. The BAO signature observed in the galaxy 
power spectrum is contaminated by the nonlinear effects 
of the density perturbations, redshift-space distortions 
and the clustering bias. The effects of the nonlinearity 
and the redshift-space distortions cause the damping of 
the BAO signature compared with the case when these 
effects are switched off. The authors of Ref. [l(| found 
that the damping of the BAO depends on the amplitude 
of the matter power spectrum and that a measurement 
of the BAO damping might be useful to constrain the 
growth factor of the cosmic structure. 

In Ref. two of the authors of the present paper 
reported a detection of the BAO damping using the Sloan 
digital sky survey (SDSS) luminous red galaxy (LRG) 
sample of the data release (DR) 6. Since the observed 
galaxy power spectrum is contaminated by errors, and 
is noisy, the procedure of extracting the BAO signature 
from a galaxy power spectrum is a subtle problem. In the 
previous paper [13], the cubic spline method was used 
to construct a smooth power spectrum to extract the 
BAO signature. However, this method is very sensitive to 
parameters of the cubic spline, i.e., number and interval 
of the nodes. Due to those reasons, it is clear that an 
independent cross check of the previous result is useful. 

In the present paper, we adopt an independent method 
to construct the smooth power spectrum, which is sta- 
ble compared with the cubic spline method, and re- 
investigate the BAO damping by confronting the model 
predictions with the SDSS LRG power spectrum of the 
DR7. We confirm the detection of the BAO damping at 
the one sigma level. Also we obtain a constraint on the 
parameter combination D\{z)<7% at the redshift z = 0.3 
from the BAO damping for the first time. Here D\(z) 
is the growth factor and erg is the amplitude of the root 
mean square overdensity in a sphere of radius 8/i Mpc. 
We denote the growth rate by f(z) = din D\(z) / d\na(z), 
where a(z) is the scale factor. Throughout this paper, we 
use units in which the velocity of light equals 1 and the 
Hubble parameter Hq = 100/ikm/s/Mpc with h — 0.7. 
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II. BAO DAMPING 

The BAO signature in the power spectrum P(fc) can 
be obtained by 



P (k) 



(i) 



where P(k) represents the 'smooth' power spectrum cor- 
responding to P(k). Since the definition of P(k) is 
not unique, the construction of P(k) should be care fully 
done. The cubic spline method is frequently used [9l.ll7j. 
but it is not always stable if dealing with a noisy power 
spectrum from observations. 

In the present paper, we construct the smooth power 
spectrum as the linear matter power spectrum multiplied 
by a simple analytic function of the wave number. With 
the use of the constructed smooth power spectrum, we 
obtain the BAO signature from the SDSS LRG power 
spectrum of the DR7 [IjJ, [2(| ■ The extracted BAO signa- 
ture is then compared with theoretical predictions, par- 
ticularly focusing on the BAO damping. 

We use the SDSS LRG sample from DR 7 (see also 
for recent results on LRGs from the SDSS DR7). 
Our LRG sample is restricted to the redshift range z = 
0.16 — 0.47. In order to reduce the sidelobes of the survey 
window we remove some noncontiguous parts of the sam- 
ple (e.g. three southern slices), which leads us to ~ 7150 
deg 2 sky coverage with a total of 100157 LRGs. The 
data reduction procedure is the same as that described in 
Ref. [6J]. In this power spectrum analysis, we adopted the 
spatially flat Lambda cold dark matter (ACDM) model 
distance- redshift relation ,s = s[z]. Figure Q] shows the 
observed power spectrum P bs(k). 
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FIG. 1: Comparison of the SDSS LRG power spectrum and 
the theoretical power spectra. The solid curve is the linear 
theory Py ln (k). The dashed curve is a 'smooth' power spec- 
trum Pobs(fc) obtained with with the Case 3. The cosmo- 
logical parameters are described in Table U labeled as model 
no.l. 
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FIG. 2: The BAO signature B(k). The solid curve is the BAO 
in linear theory B\ lrL (k). The dashed curve is the best fit the- 
oretical curve including the damping -Biin (fc) (1 — W(k)). The 
points with error bars are B b s (ki), the observational result 
obtained with the Case 3 for A(k). We adopted the cosmo- 
logical density parameters and the spectral index descried in 
Tabl^Jas model no.l. 



The observed power spectrum is contaminated by var- 
ious effects: the nonlinear effect, redshift-space distor- 
tions, clustering bias, and so on. In particular, the mod- 
eling of the clustering bias seems to be very difficult, 
because it depends on the galaxy formation process. In 
the present paper, in order to construct the smooth power 
spectrum corresponding to the observed power spectrum, 
we assume 



-fobs 

(fc) = A 2 (k)Py la (k), 



(2) 



where Pii n (fc) is the linear no- wiggle power spectrum of 
[l[ . The solid curve and the dashed curve in Fig. Q] ex- 
emplify Pnn(fc) and P bs(k), respectively. 

We determined the function A(k) as follows: We first 
define the chi-square by 

X 2 = ^[Bn„(fc l )(l-VF(fc 2 ))-Pobs(fc i )] 



1,3 



x P ohs (ki)Cov 1 (k i ,k J )P ohs (k J ) 

x [P lin (fc,)(l - W(k 3 )) - B oha {k,j)] (3) 



with 



BlinW = — — — - 1, 



-Bobs(fc) 



Pin(fc) 
Pobs(fc) 
Pobs(fc) 



(4) 
(5) 



where Cov~ 1 (fci.fc J ) is the inverse of the power spectrum 
covariance matrix [8(, and W(k) is the damping function 
(see next section for details). In the determination of 
A(k), we consider the following 3 cases, 
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FIG. 3: Ratio of the power spectrum convolved with the win- 
dow function P conv (k) to that without it P(k). The cosmo- 
logical parameters are the same as those of Fig[TJ 



Case 1 We assume A(k) = a + bk + ck 2 and determine 
a, b and c by minimising ([3]) with / = 0, <r s Di = 
(corresponding to no-damping model). 

Case 2 Same as Case 1 but with / = 0.64, a^D^ — 
0.53 (corresponding to the ACDM model with flo = 
0.28). 

Case 3 We assume A(k) = a + bk c and determine a, b 
and c by minimizing ([3]) separately for each of the 
adopted cosmologies. 

Fig [2] shows an example of the observed BAO signature 
and theoretical curves, which will be explained in the 
next section. Here the cosmological parameters are the 
same as those of Fig. [1] 

Here we note the effect of the window function, which 
we take into account in our investigation. This is the 
finite volume effect of the region where galaxies are dis- 
tributed, depending on the data analysis as well. The 
effect of the window function affects the shape of the 
BAO signature, which is not negligible in our investiga- 
tion. The details are described in the literature @, [23| . 
Fig|3] shows the ratio of the power spectrum convolved 
with the window function to that without it. 



III. THEORETICAL MODELING OF THE BAO 
DAMPING 

In our theoretical modeling, we adopt the quasi- 
nonlinear power spectrum using the technique of resum- 
ing infinite series of higher order perturbations on the ba- 
sis of the Lagrangian perturbation theory (LPT), which 
was proposed by Matsubara [H[ . One of the advantages 
of using this LPT formalism is that the redshift-space 
distortions can be taken into account. Furthermore, the 
predicted BAO signature was compared with the results 



of N-body simulations in [17j , where a good agreement 
was obtained. 

We denote the matter power spectrum in the redshift- 
space by P(k,fi), where /i is the cosine of the angle be- 
tween the line of sight direction and the wave number 
vector. However, the observed power spectrum in the 
previous section is the angular averaged power spectrum, 
which is related to P(k, n) via 



I f +1 

P(k) = - J P(k, 



fi)dfi. 



Then, the BAO signature is written as 
I-l P(k,li)dfj, 



B(k) 



(6) 



(7) 



The BAO damping is described by the function W(k), 
which is defined by 



B(k) = (1 - W(k))B nn (k) 



(8) 



where Bu n (k) is the BAO in the linear theory defined 
by Eq. (|4]). Note that B\i n (k) does not depend on the 
redshift z because the growth factor in the linear power 
spectrum cancels. We may define the angular dependent 
damping factor W(k,^L) by 



P(k,n) 
P(k,[i) 



1 



(l-IF(fc,/i))5 lin (fc), 



which leads to 



W(k) 



f 1 _ 1 dfiW(k,v)P(k,y) 
j_ 1 d^P{k,n) 



(9) 



(10) 



Following the LPT framework the matter power 
spectrum in the redshift-space is a function of /i and the 
redshift z, which we denote by -Plpt^, Mj z )- To calculate 
PhPT{k, /J,, z) one needs the linear matter power spectrum 
in redshift space, which depends on asDi(z), the growth 
factor multiplied by the normalization factor, and the 
growth rate f(z) = dlnDi(z)/dhia(z) 1 where a{z) is the 
scale factor. The full expression for -Plpt^iMi z), which 
is rather complicated, can be found in Ref. [5j. 

In the previous work [17], it was found that the an- 
gular dependent damping function at the redshift z is 
approximately given by 



W(k,n;z) = 



where Pg>(k,n) = (l + /M 2 )^F lin (fc), a(p,z) = l + /(/ + 
2)^ 2 , g(k) = (/c 2 /6tt 2 ) x J™dqPu n (q), and P {s)2 2{k^) 
is the function to describe the higher order corrections 
in the LPT framework, which is given in Appendix B of 
Ref. [5j, but note that the 'tilde' means the quantity with 
the no-wiggle linear power spectrum. 



DKz) 



(11) 
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In the present paper, we adopt W(k, fi — 0.65, z = 0.3) 
as the theoretical damping function W(k) to be com- 
pared with observations. Here z = 0.3 is a mean redshift 
of the LRG sample. Also, the angular averaged damp- 
ing function defined by Eq. (|10[) is well approximated by 
W(fc, /i = 0.65, z — 0.3), which is demonstrated in Figure 
SJ The solid curve is the angular averaged damping func- 
tion Eq. dTUD, the dotted curve is W(k, fi = 0.65, z = 0.3). 

In summary, the damping function W(k) depends on 
the linear matter power spectrum and the growth rate 
f(z) at the mean redshift z = 0.3. Then, by comparing 
the BAO damping with observations, we can hope to 
constrain a&D\(z) and f(z). 
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FIG. 4: The damping function W(fc, fi; z) at redshift z=0.3 
as a function of the wave number, for fj, = 0.5 (dot dashed), 
0.65 (dotted) and 1.0 (long dashed), respectively. The solid 
curve is the angular averaged damping function Eq. (|10|l . The 
cosmological parameters are described in Tabl^U labeled as 
model no. 1. 
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FIG. 5: The contour of A\ 2 =1.0, 2.3 on the agD^f plane 
with Case 3. The other cosmological parameters are fixed as 
those of model no. 1. in Table [I] 
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FIG. 6: The likelihood as a function of a&D\ with / integrated 
over the range < / < 2. The cosmological parameters are 
fixed as those of model no. 4. in Table [I] 



IV. COMPARISON BETWEEN OBSERVATION 
AND THEORY 

We compare the observational result and the theory 
by calculating the chi-square defined by Eq. © with the 
function A(k) determined by the methods in section 2. 
In our analysis, we treat / and a$Di as independent 
variables which have nothing to do with the cosmological 
parameters. Figure^is the contour of A\ 2 at la level on 
the a$Di-f plane, where the other parameters are fixed 
as Oo = 0.28, fib = 0.046 and the spectral index n s — 
0.96. The best fit value is a$Dt = 0.40 and / = 0.76, 
but the BAO damping is not very sensitive to the growth 
rate /. 

Figure [5] shows the likelihood as a function of og-Di, 
where the parameter / is integrated over the range < 
/<2. 



Table|T]shows the results of the chi-squared test Eq. ([3]) 
for various cosmological parameters. For the consistency 
of the analysis, we used the data of P bs(k), which was 
calculated adopting the distance redshift relation s = s[z] 
of the same cosmological parameters. From this table, it 
is clear that finite positive values of crg-Di > better fit 
the data for all the models. This allows us to conclude 
that the detection of the BAO damping is at the 1 a level 
at the worst case. The ACDM model with fJ = 0.28 
predicts D\(z = 0.3) ~ 0.65. If we adopt the value of 
er 8 = 0.8 from Ref. 0, we have a s D 1 ~ 0.53. This is 
consistent with our result within the la level. 

To discuss the attainable constraints on a$D\(z) in fu- 
ture observations, we calculate the relevant Fisher matrix 
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component 

F 1 [ kmaX d!±*( dB{k > Z) Y P ^ {k ' Z) (12) 

where we set fc max = 0.2, fc min = 0.02, B(fc, z) = (1 — 
W(fc, = 0.65), z)i?i in (fc), P ga i(/c,z) is the galaxy power 
spectrum and Q{k,z) = AA j z Zmax dz(ds/dz)s 2 n 2 /[l + 

nPg a i(fc, z)] 2 with the mean number density of galaxies n 
and the survey area AA. We adopted the Q-model for 
the galaxy power spectrum, elaborated in 12511 with the 
parameters A\ — 1.4 and Q = 16 (see also jig). 
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FIG. 7: The 1<t- level statistical errors of asDi(z) as a function 
of n(6/2) 2 with the constant number density n and the bias 
b. 



The minimum error attainable on a%D\ is given by 
l/F 1 / 2 . Figure shows 1/F 1 / 2 as a function of n{b/2) 2 , 



where & is the bias parameter of the Q-model. For the 
future surveys, we adopted the parameters in Table II jTBI, 
I2H [13] ■ We computed the constraint at the mean redshift 
z = (zmax + 2 m in)/2 for each survey. The minimum error 
CTs-Di is typically 0.1, but depends on the mean number 
density and the clustering bias. 



V. SUMMARY AND CONCLUSIONS 

In this work, we examined the BAO damping. We used 
the observed power spectrum of the SDSS LRG sample 
of DR7. The result shows that the BAO damping really 
exists and a constraint on a%Di is obtained for the first 
time from the damping. The detection of the BAO damp- 
ing is robust, and do not depend much on the specific 
model used to define the smooth model. The constraint 
on a%Di is note very tight, but this is an independent 
and unique test for the growth factor. Furthermore, this 
method might be useful for future surveys (l7l |. 

In the present analysis, we have not considered the 
effect of the clustering bias on the BAO damping. This 
point might be considered more carefully in the future. 
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No. 


iT2 m 


fit, 


n s 


f 


a$Di(la level) Case 1 


Case 2 


Case 3 


1 


0.28 


0.046 


0.96 


0.64 


0.38lg;S 


0.42l°; 3 ° 


0.421°;!* 


2 


0.28 


0.046 


0.96 


0.60 


_, nn i n 97 
U.OO_q_25 


U -^ Z -0.24 


0.42l°; 33 


3 


0.28 


0.046 


0.96 


0.70 


0.37i°;^ 


o.4olS:22 


u -^ ± -0.26 


4 


0.28 


0.046 


0.96 


integrated 


0.48lg; 3 3 ° 


_ ..-.in ^f; 
u - oz -0.30 


0.38l°; 3 * 


5 


0.27 


0.046 


0.96 


0.64 






_ _ „_i_o 4n 

0-56±g;|2 


6 


0.29 


0.046 


0.96 


0.64 


_ _ in ^9 

o.36±8:|i 


0.42l°; 3 | 


0.50l°;i 


7 


0.28 


0.048 


0.96 


0.64 


0.38lg;^ 


0-42±g;lS 


0.42l°;g 


8 


0.28 


0.044 


0.96 


0.64 


0.36lg;» 


0.42lg;g 


n 4n+°- 36 


9 


0.28 


0.046 


0.98 


0.64 


0.401°,;?,° 


0.4fit°;» 


0.421 ,; 33 


10 


0.28 


0.046 


0.94 


0.64 


0.34l°; 2 2 | 


0.38l°:S 


0.401°;^ 



TABLE I: The results of the chi-squared test of the BAO damping for various cosmological models. In the table, 'integrated' 
means that / is integrated over the range < / < 2. 



Survey 


£min 


2m ax 


n[ft 3 Mpc~ 3 ](Tyipical value) 


Area AA[deg. 2 ] 


BOSS 


0.1 


0.7 


3 x 10~ 4 


10 4 


SuMIRE 


0.7 


1.6 


4 x 10~ 4 


2 x 10 3 


EUCLID 


0.1 


2 


5 x 10~ 3 


2 x 10 4 



TABLE II: Future survey parameters adopted in the Fisher matrix analysis. 



